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1 Introduction 



Gluon fusion is the main production channel for the Standard Model Higgs boson at hadron 
colliders. Unsurprisingly, radiative corrections have been thoroughly investigated in the past years; 
in particular, since next-to-leading order (NLO) QCD corrections increase the inclusive cross section 
for Higgs production at the LHC by a factor of about 1.5 to 1.7 with respect to the leading order 
(LO) term [1], there was a flurry of activity on higher order QCD effects. Recent reviews on the 
subject can be found in Ref. [2]. 

Electroweak effects are less understood than QCD ones. An early study [3] concluded that the 
so-called leading NLO corrections to the partonic gluon- fusion cross section (j{gg — > H), enhanced 
by the squared mass of the top quark, amount to 0.4%. The size of this result is due to the 
delicate cancellation mechanism among various contributions described in Ref. [3] , and more recent 
computations have shown that these corrections are not the dominant ones in the entire Higgs mass 
range. 

Contributions induced by light quarks have been calculated in Ref. [4] and found to be extremely 
larger, reaching a maximum of about 9% for a Higgs mass below 160 GeV. The terms of the 
amplitude involving the top quark have been evaluated in Ref. [5] by means of a Taylor expansion 
in the kinematic region below the WW threshold, where they partially screen the dominant effect 
of the light quarks. The impact of the partonic results of Refs. [4, 5] on the total cross section 
for Higgs production in proton-proton collisions, a{pp H + X), has been estimated for LHC 
energies in Ref. [6]. Assuming a complete factorization of the electroweak effects with respect to 
the dominant soft and collinear QCD radiation, the analysis shows that the cross section increases 
in a range from 4 to 8 % for a Higgs mass < 160 GeV. 

A satisfactory understanding of NLO electroweak effects on Higgs production at hadron colliders 
requires at least two additional steps: first, a detailed analysis of the WW, ZZ and tt threshold 
regions, looking ahead to the possible occurrence of artificially large effects associated with the 
opening of two-particle thresholds; second, an extension of the corrections due to the top quark 
to the entire Higgs mass range. Furthermore, because of the advent of the LHC, an independent 
derivation of the results of Refs. [4-6] for a light Higgs is certainly justified. 

Improving the methods employed in Ref. [7] in the context of the Standard Model Higgs decay 
to two photons, we have recently completed the evaluation of all NLO electroweak corrections to the 
gluon- fusion Higgs production cross section at the partonic level [8] , deconvoluted of the well-known 
QCD effects [1]. 

We have found that the corrections enhance the production mechanism throughout a Higgs 
mass range spanning from 100 GeV to about 180 GeV, where light quarks dominate and the full 
NLO contributions to the partonic cross section a{gg — > H) reach up to 6%. For higher values of 
Mfj, instead, the corrections become negative and light quarks are not dominating; a minimum of 
—4 % is reached around the tt threshold. 

In addition, we have performed a dedicated study of the behavior around the WW, ZZ and tt 
thresholds [9] , showing that unphysical singularities and large threshold effects disappear once the 
complex-mass scheme of Ref. [10] is applied in a two-loop context following the strategy described 
in Ref. [11]. 

In this paper we present our numerical results for the inclusive Higgs production cross section in 
hadronic collisions, including our own evaluation of the NLO electroweak corrections in addition to 
the next-to-next-to-leading order (NNLO) QCD effects [12-14]. Moreover, we provide an estimate 
of the residual theoretical uncertainty and perform a comparison with the impact of soft-gluon 
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resummation at next-to-next-to-leading logarithmic (NNLL) level [15]. 

Our detailed numerical study is motivated by the observation that a typical size for the NLO 
electroweak corrections at the partonic level is 5 % for = 120 GeV; this value matches the 
theoretical uncertainty associated with uncalculated higher order QCD corrections, estimated to 
be 5% at the LHC and 7% at the Tevatron [16]. 



2 Inclusion of the NLO electroweak corrections 



The inclusive cross section for the production of the Standard Model Higgs boson in hadronic 
collisions can be written as 

o"(s,M^) / dxi I dx2 fi/hj {xi,f4') fj/h2 {^2,(4^) X 

„• -'0 Jo 

^ ( \ 

dz6iz ^ 1 z G,, (z; as (^|) , /^|; /4) , 
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where ^/s is the center-of-mass energy and fip and fi^ stand for factorization and renormalization 
scales. 

In Eq.([T]) the partonic cross section for the sub-process ij ^ H + X , with = g,qf,qf, has 
been convoluted with the parton densities fa/h^ ^01 the colliding hadrons hi and /12. The Born 
factor a^^^ reads 
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where Gp is the Fermi-coupling constant, Tg = M^/(4Mg) and Mq is the conventional on-shell 
mass of the top or bottom quark; the function / is 
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The coefficient functions Gij can be computed in QCD through a perturbative expansion in the 
strong-coupling constant as, 



Gij {z- as{A), Ml/fil; /^|) = a|(/x|) ^ 
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with a scale-independent LO contribution given by G 



(0), 



6ig 6jg 6{l-z). The NLO QCD 



coefficients have been computed in Ref. [1], keeping the exact Mt and Mb dependence. NNLO 
results have been derived in Ref. [12] in the large Mt limit (see Ref. [17] for the NLO case); 
analytical expressions can be found in Ref. [13] (an independent cross-check has been reported 
in Ref. [14]). The accuracy of these fixed-order computations has been improved with soft-gluon 
resummed calculations [15,16,18]. 

The inclusion of higher order electroweak corrections in Eq. ([T]) requires to define a factorization 
scheme (relevant examples on non-factorizable effects concerning Z boson decay can be found in 
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Refs. [19,20]). The authors of Ref. [6] assume that the modifications induced by sub-leading higher 
order terms starting at three loops are small; consistently with this assumption, they completely 
factorize QCD and electroweak corrections at the partonic level. This ansatz is certainly well 
justified for <^ -^w' where the electroweak interaction is effectively point-like. 

Since our analysis spans the entire Higgs mass range, and we do not foresee the advent of a 
three-loop calculation involving mixed QCD and electroweak effects, we adopt a more conservative 
approach, and resort to the well-established LEP practice of attributing a theoretical error to the 
inclusion of NLO electroweak corrections [21]. 

We introduce two options for replacing the purely QCD-corrected partonic cross section in 
Eq.([T]) with the expression including NLO electroweak corrections: 

- CF (Complete Factorization): 

^a(o) (1 + 5ew) Gij; (5) 



PF (Partial Factorization): 



G,,+a|(/z|)<^EwGSf , (6) 



where (^ew embeds all NLO electroweak corrections to the partonic cross section a{gg H), 

cTEw = a|(/ii)^^°^(l + <^Ew), (7) 

with a^^^ defined in Eq.([2]). The CF option of Eq.([5]) amounts to an overall re-scaling of the QCD 
result, dressed at all orders with the NLO electroweak correction factor S^^; the PF option of 
Eq.® is equivalent to add electroweak corrections to QCD ones. 

An intermediate option would be to fold the NLO electroweak corrections with the pure gluon- 
gluon NLO and NNLO QCD components; in this case S^w would be convoluted with Ggg^ , but not 
with c'-gj and G^^. 

Next, we define the uncertainty due to uncalculated higher order QCD corrections according 
to the standard method used in Ref. [15]: we vary the renormalization and factorization scales 
fiji and around the natural scale of the process M^j, changing their values first simultaneously, 
with the bound < fi^^p < 2M^, and then independently, with the additional constraint 

AiH/2 < !^ 2 at fixed /z^. For each value of M^, the minimal (maximal) values of the 
inclusive cross section of Eq.([T]), associated with fin,min and iJ-F,mm (A*ij,max and //F,max) are denoted 

"^min^ i^m^)', their difference defines the QCD uncertainty band around the reference value 
<T^^^ , obtained setting fi^ = fj,p = . 

Finally, for each set {M^; /i^^min; ^F.min}, {M^;na = M^;np = M^} and {M^:^ ; /i^j^max; /^F.max} 
we repeat the computation of Eq.([T]) performing the replacements of Eq.(l5|) and Eq.Q. As a 
result, two new reference values are obtained for both CF and PF options, giving the impact of our 
evaluation of the NLO electroweak corrections. In addition, we obtain new minimal and maximal 
values for the cross section, dressed with electroweak effects, and define a new uncertainty band. 

Needless to say, this procedure will only give an approximate bound on the total uncertainty; 
one should indeed observe that for the considered range the LO and NLO QCD bands of Ref. [1] 
do not overlap, with an NNLO band which is only partly contained in the NLO one. 
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3 Numerical results 



For the NLO electroweak corrections we use our recent result [8] and consider a Higgs mass 
range spanning from 100 GeV to 500 GeV. In this region we cross the WW, ZZ and thresholds. 
A naive computation of the amplitude with conventional on-shell masses as input data reveals the 
presence of singularities at the WW and ZZ thresholds; in order to cure them, we have introduced 
in our computation complex masses [11], following the suggestion of Ref. [10]. The behavior at the 
ti thresholds, instead, is smooth, and the on-shell mass of the top quark can be safely used. 

In the calculation all light-fermion masses have been set to zero and we have defined the W and 
Z boson complex poles by 

Sj = fij {fij - i 7j) , fi] = M] - V], 7j = Tj ^1 - , (8) 

with j = W, Z. As input parameters for the numerical evaluation we have used the following values 
taken from Ref. [22]: 

= 80.398 GeV, = 91.1876 GeV, , 

= 2.4952 GeV, = 1.16637 x 10"^ GeV^^. 

For the mass of the top quark, we have used Mt = 170.9 GeV [23]; for the width of the W boson, 
we have chosen the value = 2.093 GeV, predicted by the Standard Model with electroweak and 
QCD corrections at one loop. 

Our results for 5ew defined in Eq.© are shown in Fig. [H where we include the complete 
corrections, comprehensive of light- and top-quark contributions, in the entire range of interest. The 
introduction of the complex-mass scheme in our two-loop evaluation has a striking consequence, 
visible around two-particle thresholds, where artificial cusp effects disappear. A detailed analysis 
of this issue can be found in Ref. [9] . 

6 
4 

^ 
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Figure 1: NLO electroweak percentage corrections to the partonic cross section a{g g H). 

For including the NLO electroweak corrections in the hadronic cross section of Eq.([T]), we have 
used the FORTRAN code HIGGSNNLO written by M. Grazzini (see also Ref. [24]), with QCD 
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corrections at NNLO, interfaced with the MRST2002 set of parton distribution functions [25]. 
Although partially outdated, they represent the best choice for our purposes, allowing for a direct 
comparison with the results of Ref. [15]. In the code the parton densities, as well as the strong- 
coupling constant Os, are evaluated at each corresponding order, with one-loop as at LO (as(M^) = 
0.130), two-loop as at NLO (as(M2) = 0.1197) and three-loop as at NNLO {asiM^) = 0.1154), 
as described in Ref. [15]. 

In the following, we will discuss our numerical results introducing K factors, defined as the 
ratio of the cross section including higher order corrections over the LO result. According to the 
discussion in Section [21 we will define K factors for NNLO QCD corrections and for NNLO QCD 
-|- NLO electroweak corrections under the assumption of complete (partial) factorization. The LO 
expression which normalizes each K factor will be always evaluated for fi^ = fip = , with parton 
densities and as evolved at LO. 

3.1 LHC results 

We start showing our results for a LHC center-of-mass energy ^/s = 14 TeV. According to 
the standard procedure described in Section [21 an uncertainty band for the K factors is derived 
exploring the dependence of the cross section on the renormalization and factorization scales. In 
Fig. [21 we show the maximal and minimal values as a function of , for 100 GeV < < 500 
GeV, obtained including NNLO QCD corrections only (dashed hues) and NNLO QCD + NLO 
electroweak ones (solid lines). Note that in the second case the maximal and minimal values take 
into account both factorization options of Eq.(l5]) and Eq.([6]). The pattern of the electroweak- 
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Figure 2: Uncertainty bands for the K factors for Higgs production at the LHC. 



pp ^ H + X Vs = 


14 TeV ^^.^is-a^ 


MRST 2002 




/ \ ^."^ 










— — ■ — 1 y^ 







/- — 




/ y 


NNLO QCD 

NNLO OCD + NLO EW 



corrected uncertainty band for the K factors in Fig. [2] clearly follows the shape of the correction 
factor at the partonic level of Fig. [H At about 180 GeV, when (5ew becomes negative, the two 
factorization options exchange their role: below 180 GeV, the completely (partially) factorized 
result fixes the upper (lower) bound, above 180 GeV the completely (partially) factorized result 
fixes the lower (upper) bound. This is just a consequence of the fact that the assumption of a 
complete factorization leads to a larger enhancement of the absolute value of the NNLO QCD 
corrections. 
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In contrast with the purely QCD-corrected K factors, we observe a higher sensitivity of the 
electroweak-corrected shapes of Fig. [2] on the values of the Higgs mass. The larger effect takes 
place at = 150 GeV, where the QCD band at NNLO ranges from 1.92 to 2.33 and NLO 
electroweak effects shift it to an interval between 1.97 and 2.47. 

The corresponding numerical values for the NNLO QCD-corrected cross section a^*-'^, and 
the cross section including NNLO QCD + NLO electroweak corrections under the assumption of 
complete (partial) factorization, {a^^)^ are shown in Tab. [H including minimal, reference and 
maximal values (the differences with Ref. [15] are simply due to the different values used for the 
top-quark mass, Mt = 170.9 GeV in this paper, Mt = 176 GeV in Ref. [15]). 

Strictly speaking, the values for the QCD corrections at NNLO are given here in the large- Mt 
limit; however, it has been shown at NLO that this approximation is extremely good for < 2 Mt 
and works with a good accuracy up to M^^ = 1 TeV [26]. 
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100 


47.84 


53.44 


59.18 


49.83 


55.68 


61.65 


48.73 


54.55 


60.57 


110 


41.22 


45.92 


50.72 


43.09 


48.01 


53.02 


42.05 


46.94 


51.99 


120 


35.94 


39.96 


44.03 


37.71 


41.92 


46.19 


36.72 


40.91 


45.21 


130 


31.66 


35.13 


38.63 


33.34 


36.99 


40.68 


32.40 


36.02 


39.74 


140 


28.14 


31.16 


34.20 


29.73 


32.92 


36.14 


28.83 


31.99 


35.24 


150 


25.20 


27.86 


30.53 


26.69 


29.50 


32.33 


25.84 


28.63 


31.49 


160 


22.73 


25.08 


27.45 


23.82 


26.29 


28.77 


23.19 


25.64 


28.15 


170 


20.62 


22.73 


24.84 


21.03 


23.18 


25.33 


20.79 


22.94 


25.10 


180 


18.82 


20.72 


22.61 


18.91 


20.81 


22.72 


18.86 


20.76 


22.67 


190 


17.27 


18.99 


20.70 


17.04 


18.74 


20.43 


17.18 


18.87 


20.56 


200 


15.93 


17.49 


19.05 


15.59 


17.12 


18.65 


15.79 


17.32 


18.85 


220 


13.75 


15.07 


16.39 


13.41 


14.70 


15.98 


13.61 


14.91 


16.18 


240 


12.11 


13.25 


14.38 


11.81 


12.93 


14.03 


11.99 


13.10 


14.21 


260 


10.87 


11.87 


12.87 


10.63 


11.61 


12.58 


10.77 


11.76 


12.73 


280 


9.96 


10.87 


11.76 


9.76 


10.65 


11.53 


9.88 


10.77 


11.65 


310 


9.15 


9.98 


10.80 


8.99 


9.81 


10.61 


9.09 


9.91 


10.70 


340 


9.38 


10.24 


11.07 


9.15 


9.98 


10.79 


9.29 


10.13 


10.93 


370 


10.50 


11.46 


12.39 


10.16 


11.10 


11.99 


10.37 


11.31 


12.19 


410 


9.00 


9.83 


10.62 


8.83 


9.65 


10.42 


8.94 


9.75 


10.52 


450 


6.91 


7.55 


8.15 


6.86 


7.49 


8.09 


6.89 


7.52 


8.12 


500 


4.73 


5.17 


5.58 


4.74 


5.18 


5.59 


4.73 


5.17 


5.59 



Table 1: NNLO QCD and NNLO QCD + NLO electroweak (with CF (Eq.©) and PF (Eq.®) 
options) cross sections in pb as a function of the Higgs mass in GeV at the LHC. 
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3.2 Tevatron results 

In this section we briefly summarize the results for the Tevatron Run II, with a center-of- 
mass energy -y/s = 1.96 TeV. In Fig. [3] we show the maximal and minimal values for the K 
factors as a function of in the range 100 GeV < < 200 GeV, obtained including NNLO 
QCD corrections (dashed lines) and NNLO QCD + NLO electroweak ones (solid lines). As usual, 
maximal and minimal values take into account both factorization options of Eq.([5]) and Eq.([6]). 
The corresponding numerical values are shown in Tab. [2j 

As stressed in Ref. [15], the K factors at the Tevatron are larger than those at the LHC, 
because the Higgs production takes place closer to the hadronic threshold and soft-gluon effects 
are extremely relevant. Consequently, the impact of soft-gluon resummation is more sizeable than 
the electroweak effects: for = 120 GeV, the NNLL result of Ref. [15] increases the fixed-order 
NNLO value by 12 %; from Tab. [H the impact of NLO electroweak terms amounts to +5 % (+2 %) 
assuming a complete (partial) factorization. 

Note that the same consideration is not true for the LHC, where NLO electroweak effects are 
comparable to those due to soft-gluon resummation at NNLL accuracy. For = 120 GeV, we 
get from Ref. [15] an increase of 6 %, which matches the electroweak effects; the latter amount to 
+5% (+2%) assuming a complete (partial) factorization (see Tab.[l]). Further details will be given 
in Section 13.41 
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Figure 3: Uncertainty bands for the K factors for Higgs production at the Tevatron. 



3.3 Threshold behavior 

The scheme employed in our computation [8,9] shows that the partonic cross section (j{gg — > H) 
is a smooth function of the Higgs mass; as illustrated in Fig. [H there are no artifical large effects 
at the opening of two-particle thresholds. 

As a consequence, unphysical cusps are avoided also at the hadronic level; in Tab. [3] we sum- 
marize the results at the LHC for values of the Higgs mass corresponding to the WW , ZZ and ti 
thresholds, following the same pattern of Tab. [T]and Tab. [2j 
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QCD 
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QCD 
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rcf 
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100 


1.2981 


1.4886 


1.6829 


1.3523 


1.5508 


1.7532 


1.3136 


1.5097 


1.7122 


110 


1.0002 


1.1482 


1.2978 


1.0455 


1.2003 


1.3567 


1.0130 


1.1656 


1.3222 


120 


0.7833 


0.9002 


1.0175 


0.8219 


0.9445 


1.0675 


0.7940 


0.9147 


1.0380 


130 


0.6218 


0.7153 


0.8086 


0.6548 


0.7533 


0.8515 


0.6307 


0.7276 


0.8260 


140 


0.4993 


0.5750 


0.6500 


0.5276 


0.6075 


0.6868 


0.5068 


0.5854 


0.6648 


150 


0.4049 


0.4668 


0.5279 


0.4288 


0.4943 


0.5590 


0.4112 


0.4755 


0.5403 


160 


0.3314 


0.3824 


0.4326 


0.3473 


0.4008 


0.4534 


0.3355 


0.3882 


0.4408 


170 


0.2734 


0.3158 


0.3573 


0.2788 


0.3220 


0.3644 


0.2748 


0.3177 


0.3601 


180 


0.2272 


0.2627 


0.2974 


0.2283 


0.2639 


0.2988 


0.2275 


0.2630 


0.2979 


190 


0.1901 


0.2200 


0.2492 


0.1876 


0.2171 


0.2459 


0.1895 


0.2191 


0.2479 


200 


0.1601 


0.1854 


0.2101 


0.1567 


0.1815 


0.2057 


0.1593 


0.1843 


0.2084 



Table 2: NNLO QCD and NNLO QCD + NLO electroweak (with CF (Eq.©) and PF (Eq.®) 
options) cross sections in pb as a function of the Higgs mass in GeV at the Tevatron. 
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160.8 


22.55 


24.88 


27.23 


23.58 


26.02 


28.47 


22.99 


25.41 


27.88 


182.3 


18.44 


20.29 


22.14 


18.44 


20.29 


22.14 


18.44 


20.29 


22.14 


341.8 


9.50 


10.37 


11.21 


9.21 


10.05 


10.86 


9.39 


10.23 


11.04 



Table 3: NNLO QCD and NNLO QCD + NLO electroweak (with CF (Eq.©) and PF (Eq.®) 
options) cross sections in pb at the two-particle thresholds at the LHC. is given in GeV. 

At the WW threshold, electroweak effects amount to +5% (complete factorization, CF) and 
+2% (partial factorization, PF) of the NNLO QCD result; around the ZZ threshold, they are 
vanishingly small, since here the corrections to the partonic cross section of Fig. [1] are negligible. 
Finally, the effect at the ti threshold is to moderately decrease the cross section, by an amount of 
-3% (CF) and -1% (PF). 

3.4 Comparison with NNLL soft-gluon resummation 

In Tab. H] we show a comparison, for ^/s = 14 TeV, of our result with the NNLL resummation 
performed by the authors of Ref. [15]. For the NNLO QCD cross section (with and without NLO 
electroweak corrections) we define an average value, cxaver = ('7max+Cmin)/2, and the associated error 
Ao" = ((Tmax — Cmin)/2. Maximal and minimal values have been shown in Tab. [H for electroweak 
corrections, we have taken into account both CF and PF options. The resulting electroweak shift 
can be directly derived from the second and third columns of Tab. HI the NNLL shift is obtained 
from Tab. 1 of Ref. [15], taking the difference of the two reference values for the fixed-order and 
resummed computations. 
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200 


17.49 


± 


1.56 


17.22 ± 


1.63 


- 0.27 


+ 


1.04 


220 


15.07 




1.32 


14.80 ± 


1.39 


- 0.27 


+ 


0.89 


240 


13.24 




1.14 


13.01 ± 


1.20 


- 0.23 


+ 


0.79 


260 


11.87 




1.00 


11.68 ± 


1.05 


- 0.19 


+ 


0.72 


280 


10.86 




0.90 


10.70 ± 


0.94 


- 0.16 


+ 


0.65 



Table 4: Shifts induced at the LHC on the NNLO QCD cross section by NLO electroweak effects 
and NNLL resummation taken from Ref. [15]. The Higgs mass is in GeV, all values for the cross 
section and the shifts are in pb. 

We observe that for < 160 GeV the inclusion of NLO electroweak corrections increases the 
average value of the cross section by about half the size of the QCD uncertainty half-band; the effect 
is comparable with the impact of soft-gluon resummation at NNLL accuracy. For higher values 
of the Higgs mass, NLO electroweak effects are negative, and they partially screen the impact of 
soft-gluon resummation on the fixed-order NNLO QCD result. 

4 Conclusions 

In this paper we have presented the impact of the NLO electroweak corrections to the inclusive 
cross section for Higgs production at the LHC. 

The central value of the NNLO QCD cross section for = 120 GeV is shifted by +5 %, 
both at the LHC and the Tevatron, under the assumption of a complete factorization of the NLO 
electroweak effects with respect to the dominant QCD radiation. The impact is relevant in view of 
the estimated uncertainty associated with uncalculated higher order QCD corrections, 5 % at the 
LHC and 7% at the Tevatron [16]. The underlying assumption is motivated by the observation 
that for low Higgs masses the electroweak interaction is effectively point-like [6]. We have also 
derived a more conservative estimate on the electroweak shift assuming a partial factorization at 
leading order: for = 120 GeV the result is reduced to +2 %. 

The NNLO QCD fixed-order uncertainty bands, derived varying the renormalization and factor- 
ization scales at fixed values for the Higgs mass, have been refined including the NLO electroweak 
effects; as a result, they show a stronger sensitivity to the Higgs mass with respect to the pure 
QCD result. 

For low values of the Higgs mass, we get a qualitative agreement with the results of Ref. [6], as 
a consequence of the light-quark dominance. For higher values of the Higgs mass, the role of the 
top quark becomes relevant and the agreement starts deteriorating. 
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We have shown that large two-particle threshold effects are avoided working in the complex-mass 
scheme of Ref . [10] . Further details can be found in our companion paper [9] . 

We have performed a detailed analysis for a wide Higgs mass range, up to 200 GeV at the 
Tevatron and 500 GeV at the LHC. Concerning the Tevatron, we have confirmed the expectations 
of Ref. [15], showing that electroweak effects are considerably smaller than those induced by the 
soft-gluon resummation. At the LHC, instead, the size of the positive NLO electroweak corrections 
is comparable to that of the positive soft-gluon resummation at NNLL; above 180 GeV, they are 
negative and moderately screen soft-gluon effects. 

In summary, electroweak effects to Higgs production at hadron colliders are under control at 
NLO in the whole Higgs mass range. The main source of uncertainty is connected with a more 
precise knowledge of the parton distribution functions. 
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